Entanglement and quantum phase transition in the extended Hubbard model 
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We study quantum entanglement in one-dimensional correlated fermionic system. Our results 
show, for the first time, that entanglement can be used to identify quantum phase transitions in 
fermionic systems. 
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Quantum entanglement, as one of the most intrigu- 
ing feature of quantum theory, has been a subject of 
much studies in recent years, mostly because its non- 
local connotation^ is regarded as a valuable resource 
in quantum communication and information process- 
ing H Q One important issue is whether there ex- 
ists any relation between quantum entanglement and 
quantum phase transitions 0. Several groups inves- 
tigated this problem by study quantum spin systems 
[I II II EMI El 111 El- For examples the work of 
Osterloh et aZ.@] and Osborne and Nielsen)^] on the spin 
model showed that the entanglement of two neighboring 
sites displays a sharp peak either near or at the critical 
point where quantum phase transition undergoes. 

On the other hand, real systems consist of moving 
electrons with spin so to explore the relation between 
quantum entanglement and quantum phase transition in 
fermionic system is necessary. Previously, there are cou- 
ple of works studied entanglement in fermionic lattices 
|l5l Il6| . but they did not discuss its relation to quan- 
tum phase transition. In this Letter, in the framework 
of one-dimensional extended Hubbard model, we study 
the change of symmetry in the ground state on pass- 
ing the phase boundary from the point view of quantum 
entanglement, and demonstrate that entanglement is an 
unique quantity to describe quantum phase transitions 
in this system. The one-dimensional extended Hubbard 
model (EHM) is defined by the Hamiltonian 
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In Eq. Q, a =t, |; j — 1,...,L;S = ±1, c) 
are creation and annihilation operators of electron with 
spin a at site j, respectively. U and V define the on-site 
and nearest-neighbor Coulomb interactions. The EHM 
is a prototype model in condensed matter theory for it 
exhibits a rich phase diagram 0, 0, 0] where various 
quantum phase transitions occur between symmetry bro- 
ken states. These states include the charge-density-wave 
(CDW), the spin-density- wave (SDW), and phase sep- 
aration (PS). By calculating the entanglement as func- 
tions of electron-electron interaction U and V as well 
as fermion concentration N/L, we show that quantum 
phase transitions can be identified at places where local 



entanglement is extremum or its derivative is singular. 
Our results, part of which are based on the exact solu- 
tion of the one-dimensional Hubbard model, arc useful 
for people to explore quantum entanglement and quan- 
tum phase transition via other approaches for interacting 
many-fermion systems. 

For spin-1/2 fermion system, there are four possible 
local states at each site, \u)j = \Q)j, | t)j, | \ 
The dimension of the Hilbcrt space for a L-site system 

is then A L , and v% • • • vl) = Ilj=r \ v j)i are its nat- 
ural basis vectors. We consider local density matrix of 
the ground state, |\&), which is a reduced density matrix 
Pj = Trj | >F) (>F | , where Trj stands for tracing over all sites 
except the jth site. Accordingly, the von Neumann en- 
tropy E v calculated from the reduced density matrix pj 
measures the entanglement of states on the jth site with 
that on the remaining L — 1 sites. It is called the local 
entanglement |16| for it exhibits the correlations between 
a local state and the other part of the system. Since 
Hamiltonian Q is invariant under translation, the local 
density matrix pj is site independent, 

Pj = z | o)(o | +u+ ITXTI +u~ |j)UI +w IUXUI , (2) 

with 

w = {n n n n ) = Tr(n n n jlPj ), 

u + = (n T )-u>, u" = (n i )-u;, 

z = 1 — u + — u~ — w = 1 — (rtj) — (rtj) + w . (3) 

Consequently, the corresponding von Neumann entropy 
(or the local entanglement which we call hereafter) 

E v = —z log 2 z — u + log 2 u + — u~ log 2 u~ — w log 2 w 

Clearly, the local entanglement combines four quantities 
which arc all important to decide the physical properties 
of the system. We discover, to be shownn below, that 
this simple expression plays more general and important 
role for the understanding of the system than other single 
parameter. 

We start with general behavior of the local entangle- 
ment for the half- filling (N = L) case. In Fig. 1, we 
plot E v on the U — V plane with its contour map. It is 
remarkable to see that the skeleton of the EHM's phase 
diagram |19| can be directly obtained from the contour 




FIG. 1: The changes of symmetry in the ground state wave- 
function is analyzed by considering the quantum correlation 
between local site and other part of the system. The curved 
surface denotes E v 's dependence on U and V, and colored 
curves on E v — plane constitutes a contour map. Three 
solid lines on the plane denote the the local extremum of a 
transact of "mountain" surface. Clearly, three main symme- 
try broken phases (CDW, SDW and PS) can be sketched out 
from the contour map. Superconducting phase could not be 
identified, due to the fact that the broken symmetry is asso- 
ciated with off-diagonal long-range-order. Finite size scaling 
analysis should be carried out. 



map. This is by no means trivial. In conventional ap- 
proach to obtain the phase diagram of the EHM, one has 
to study behaviors of different order parameter in differ- 
ent regions, either by comparing ground state energy or 
critical exponent of correlation function associated with 
broken symmetry. Whereas here, using a single quantity, 
E v , the global picture of the system could be observed. 
Just like the proverb says, "a drop of water reflects the 
rays of the sun". Obviously, this is not a coincident. 
Rather, it reflects the underlying correlation between en- 
tanglement and quantum phase transition behind the su- 
perposition principle of quantum mechanics. 

In order to clarify physical pictures further, we present 
our studies in details at some special transects. Firstly, 
we study the Hubbard model, i.e., V = in the EHM. 
Since the one-dimensional Hubbard model can be solved 
analytically for both finite and infinite lattices by the 
Bethe-ansatz method 201. we can study the analyticity 
of the phase transition as well as checking the validity of 
the numerical exact diagonalization technique used for 
the EHM on finite lattices. 

For the Hubbard model it is well known that its ground 
state is a spin singlet and (nj) = (nj,) = 1/2. Thus, u + 
and u~ are both equal to 1/2 — w and the local entan- 
glement is 



FIG. 2: Local entanglement E v of the Hubbard model at half- 
filling versus the on-site coupling U for different size lattice. 



By making the use of particle-hole symmetry of the 
model, one easily finds that w{—U) = 1/2 — w(U), so 
the local entanglement is an even function of U, i.e. 
E v (—U) = E V (U). In the large U limit, \U\ — > oo, either 
all sites are singly occupied (U > 0) so w — 0, or half of 
the total sites are doubly occupied while the other half 
are empty sow= 1/2, one gets E V (\U\ = oo) = 1. For 
finite \U\, hopping process enhances E v , which reaches 
its maximum value 2 at U = from both sides. We plot 
the local entanglement E v as functions of U in Fig. 
obtained from the Bethe ansatz method (for L = oo and 
L = 70) and exact diagonalization technique (L = 10). 
The excellent agreement justifies the validity of using 
small clusters for other calculations. The ground state 
of the one-dimensional Hubbard model at half-filling is 
metallic for U < 0, and insulating for U > 0, so U = 
is a critical point which separates metallic and insulating 
phases. Thus, our result shows that the local entangle- 
ment reaches its extremum at the critical point where the 
system possesses maximum SO (4) symmetry and under- 
goes quantum phase transition. 

Moreover, based on the Bethe ansatz solution, we can 
also study the asymptotic behavior of the entanglement 
analytically. In the large U 3> 1 region, to the third or- 
der in l/C/ 2 [2ll, we have w = 41n2/[/ 2 - 27((3)/U 4 + 
375C(5)/U 6 where £ stands for the Riemann zeta func- 
tion. Therefore the local entanglement yields the fol- 
lowing asymptotic behavior E v = 1 + 161n[//[/ 2 + • ■ • . 
Whereas in the week coupling region < U <C 1, 
the density of double occupancy becomes w = 1/4 — 
7C(3)(7/87r 3 -93C(5)t/ 3 /2 9 7r 5 , which is obtained bymak- 
ing the use of energy expansion with respect to U |22ll23| . 
Thus, the local entanglement near the critical point is, 



E v = 2 — 
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K(,3)U 
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E,, 



-2wlog 2 w - 2 (1/2 - w) log 2 (1/2 - w) . (4) 



Clearly, E v is analytic in 

the neighborhood of the critical point U = 0. This be- 
havior is different from other models 0, 0, ^| . 



Secondly, we study the EHM at some fixed values of on- 
site Coulomb interaction U by varying nearest-neighbor 
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FIG. 3: Local entanglement E v versus V for some values of 
U. 
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FIG. 4: Local entanglement E v as a function of fermion 
concentration for the Hubbard model with different on-site 
!7(left), and the first derivative of E v with respect to n at 
n — 1 as a function of U. 



interaction V. For negative U, each site tends to be dou- 
bly occupied. When V ^> \t\, CDW state is favored, while 
for V <C — |t|, phase separation occurs. Both CDW state 
and PS state lead to E v = 1. Only in the region where 
\V\ ~ 0, electron itinerant motion dominates, tends to 
uniform density distribution so a local maximum of E v 
occurs around V = 0, as shown in Fig. |3] 

For positive U, the physics becomes more interesting 
since the model in this region is more relevant to the real 
materials. When V > 0, the competition between CDW 
and SDW will lead E v to an extremum where the phase 
transition undergoes, due to the fact that the local entan- 
glement itself combines CDW order parameter and SDW 
order parameter at the same time. As shown in Fig. 
121 and |21 the transition happens along a line U rs 2V, 
consistent with other studies of the EHM [l7L HI flQj . 
When V < 0, the formation of electron clustering, i.e., 
the phase separated configuration, challenges the SDW 
state. In the large U and \ V\ limit, it can be easil y sh own 
that phase transition happens at U = —2V |l7L Hsl fl9|. 
^From Fig. |3 we also observe that there exists singu- 



lar behavior for the local entanglement at the transition 
points. Moreover, for U > 0, we find that the local en- 
tanglement at two boundary lines is very close to 2, indi- 
cates that each of the four local modes has nearly equal 
population at critical point. 

Thirdly, we study the variation of local entanglement 
as function of chemical potential by adding the term 
—/j, rii to the Hubbard model. Consequently, the total 
particle number of the ground state, hence the filling fac- 
tor, could be tuned. We show the relations between local 
entanglement and the filling factor n for various on-site 
coupling U in Fig. 0] We only need to plot the part of 
n = N/L < 1 because the other part, n > 1, could sim- 
ply be obtained by the mirror image relation, as easily 
seen by the particle-hole transformation, namely, 



E v (n)=E v (2-n). 



(5) 



Fig. 01 manifests that the ground state of the half-filled 
band is not maximally entangled as long as U > 0, 
whereas, the maximum of E v lies in between n = 2/3 
and n=l. Let us take U = oo for example. When 
U = oo, there is no double occupied site, which im- 
plies that w — and u + = u~ = N/2L. Hence we 
have an analytical expression of the local entanglement 
E v — —(1 — n)log 2 (l — n) — nlog 2 (n/2) which has a 
maximum at n = 2/3. It is worthwhile to point out that 
at 1/3 filling (i.e., n = 2/3) when U = oo, the ground 
state is a singlet of SU(2|1) Lie supersymmetry algebra 
which possesses the maximal symmetry allowed, while at 
1/2 filling, it is a SU(2) singlet. For U = the ground 
state is invariant under SO (4) rotation at 1/2 filling. This 
demonstrates that the local entanglement reaches a maxi- 
mum value at the state with maximal symmetry. Accord- 
ingly, the maximum position for < U < oo is expected 
to lie between n = 2/3 and n — 1, which is numerically 
confirmed in Fig. 01 

Except at half-filling where it becomes a Mott- 
insulator, the system is an ideal conductor [2(|. Conse- 
quently, the local entanglement E v is not smoothly con- 
tinuous at n = 1 for U ^ 0. It is then instructive to 
observe the derivative of E v with respect to U, 



dE v 
dn 



(log 2 u + - log 2 z) 



71 = 1- 



1 dAE 
- + 2 

2 dU 



,(6) 



71=1 



where AE is the gap of charge excitation. Eq.lJSJ gives 
rise to dE v /dn\ n=1+ = — dE v /dn\ n=1 - . Obviously, 
there exists a jump in the derivative of E v across the 
point of insulating phase (see Fig. QJright) unless U = 0. 

^,From the above investigations, we find that the lo- 
cal entanglement manifests distinct features at the point 
where quantum phase transition undergoes. Since the lo- 
cal entanglement represents the symmetry of the system 
to a certain extent, naturally one expects that the maxi- 
mum point of the local entanglement corresponds to the 
quantum phase transition point. In the light of this con- 
clusion, we speculate that the maximum point in Fig. 01 
not only denotes the maximum symmetry, but could also 
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be a critical point separating two different phases. On 
the other hand, the discontinuity properties of the local 
entanglement obviously indicates a phase transition. For 
example, the derivative of E v in the region of U > and 
V < at half-filling of the EHM, and of the Hubbard 
model caused by the shifting of the chemical potential 
are both not smoothly continuous at the quantum phase 
transition points. This is similar to other studies, e.g., 
the one-dimensional XY model in a transverse magnetic 
field iiia, where the derivative of the pairwise con- 
currence C with respect to the dimensionless coupling 
constant develops a cusp at the quantum phase transi- 
tion point. However, such discontinuity is not universal, 
as shown by our results. 

It was indicated [24J recently that two mechanisms 
may bring about quantum phase transitions in one- 
dimensional correlated fermionic systems. One is caused 
by the level crossing of the ground state and the other 
arises from the level crossing of the low-lying excited 
states where no level crossing occurs at the ground state. 
In the later case, if the ground state wavefunction is 
smoothly continuous with respect to the variation of pa- 
rameters that drive the quantum phase transition, the 
entanglement should also be smoothly continuous at the 
quantum phase transition point. On the other hand, the 
singularity of wavefunction may leads to the singularity 
of the local entanglement. For the former case, the level 
crossing of the ground state will clearly cause the entan- 
glement to be none smoothly continuous at the transition 
point. Therefore the continuity property of the local en- 
tanglement might be an ancillary tool to judge the mech- 
anism of quantum phase transition proposed in |24| . 

In summary, we have extensively studied the local en- 



tanglement in the one-dimensional extended Hubbard 
model, characterized by the on-site Coulomb interac- 
tion U, the nearest-neighbor Coulomb interaction V , and 
band filling N/L. At half filling, we calculated local en- 
tanglement as functions of U and V. Our results indi- 
cated that the local entanglement either reaches the max- 
imum value or shows singularity (or both) at the critical 
point where quantum phase transition undergoes. For 
the traditional Hubbard model (V = 0), the scaling be- 
havior close to the critical point [7 = was given as 
manifests that the local entanglement is an analytical 
function of U. The asymptotic behavior of the local en- 
tanglement at the strong coupling limit, U — ► oo, was 
also given. Furthermore, we analyzed the local entangle- 
ment by varying nearest-neighbor V while keep U fixed, 
and found that the local entanglement is not smoothly 
continuous in some critical regions. Finally, we studied 
the dependence of the local entanglement on the filling 
factor for the Hubbard model. The variation of the lo- 
cal entanglement caused by the shifting of the chemi- 
cal potential showed that the local entanglement reaches 
maximum at filling factor n between 2/3 and 1. For 
any finite U, the local entanglement develops a cusp at 
n = 1. In the strong coupling limit, the 1/3 filled band 
has the maximum local entanglement, suggesting that 
the ground state with maximal symmetry possesses the 
maximum magnitude of the local entanglement. 
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